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Abstract 



We calculate microscopically total and parity-projected level densities for 
/3-stable even-even nuclei between Fe and Ge, using the shell model Monte 
Carlo methods in the complete (pf + Ogg^-shell. A single-particle level den- 
sity parameter a and backshift parameter A are extracted by fitting the cal- 
culated densities to a backshifted Bethe formula, and their systematics are 
studied across the region. Shell effects are observed in A for nuclei with 
Z = 28 or N = 28 and in the behavior of A/ a as a function of the number 
of neutrons. We find a significant parity-dependence of the level densities for 
nuclei with A < 60, which diminishes as A increases. 
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Nuclear level densities are important in theoretical estimates of various compound nuclear 
reactions. For example, neutron-capture reaction rates are approximately proportional to 
the corresponding level densities in the neutron resonance region. Neutron capture rates are 
relevant to the s- and r-processes in nucleosynthesis. The s-process proceeds by a sequence of 
neutron-capture reactions along the /5-stable line, since the /3-decay of a radioactive nucleus 
will occur much faster than the capture of a neutron. The waiting points of the r-process 
are determined by the balance between the rates of neutron-capture and of photoejection of 
a neutron. Accurate methods for calculating level densities are especially useful for nuclei 
where neutron-capture cross-sections are unavailable experimentally. Nuclei in the iron- 
group are of special interest since they are the seed nuclei for the synthesis of the heavy 
elements by the s- and r-processes. 

Experimentally, level densities are usually extracted from direct counting at low exci- 
tation energies (E x < 5 MeV), from neutron or proton resonance data |]], from charged 
particles spectra at intermediate energies (5 < E x < 15 MeV) [@,[|, and from Ericson fluc- 
tuation analysis at higher energies (E x > 15 MeV) @. In many nuclei, the densities at 
low energies and around the neutron resonance region are well-fitted to phenomenological 
modifications of the Fermi gas model [|j], in particular to the backshifted Bethe formula 
(BBF) [||||. However, the fitted single-particle level density parameter a and backshift pa- 
rameter A are nucleus- dependent and difficult to derive theoretically. Consequently, level 
densities used in nuclear astrophysical studies are based on various empirical formulae for 
a and A |5|-f§. Recently we introduced a novel method to calculate microscopically total 
and parity-projected level densities |J in the framework of the auxiliary-field Monte Carlo 
method. The method takes into account exactly correlations due to effective two-body in- 
teractions in a finite shell model space. We showed that realistic nuclear level densities can 
be calculated with good accuracy if a sufficiently large model space is used. In particular, 
we applied the method to 56 Fe and found the calculated total level density to be in good 
agreement with the experimental one. We found a significant parity-dependence of the level 
density, contrary to the common assumption in astrophysics calculations that positive- and 
negative-parity levels have equal densities around the neutron resonance region. Shell (or 
subshell) structure affects the parity dependence of level densities, as pointed out already 
in Ref. [T(J for non-interacting fermions. In this letter, we investigate the systematics of 
the total and parity-projected level densities of /3-stable even-even nuclei between Fe and 
Ge. The calculated level densities are found to be well-described by the BBF. We extract 
both a and A from the microscopic Monte Carlo densities by fitting them to the BBF, and 
find shell effects in their systematics. The extracted parameters are also compared with 
empirical values and with experimentally extracted values when available. 

Approximate methods to calculate level densities beyond the Fermi gas model were de- 
veloped in the framework of the spectral averaging theory |~lHIl]- Level densities for several 
p/-shell nuclei were recently calculated [|13j by decomposing the level density into a sum of 
non-interacting level densities convoluted with Gaussians (the Gaussians reflect the effect of 
interactions). Since Gaussian spreading is justified only above a certain reference energy, the 
method used in Ref. |13|] requires the knowledge of a complete set of levels up to sufficiently 
high energy. While reasonable agreement with data is found, two parameters, both of which 
depend on the nucleus under study, were fitted to the level density data: a reference energy 
and an interaction parameter. In our level density calculations, no adjustable parameters 
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have been introduced, yet the calculated level densities are in good agreement with the data. 

Neutron resonance energies for 50 ^5 A < 70 nuclei are in the range E x < 15 MeV. To 
include the important excitations in this mass and energy region in a shell model approach, 
it is necessary to use the full (pf + 0gg/ 2 ) model space. This model space is too large for 
conventional diagonalization methods: instead we use the shell model Monte Carlo (SMMC) 



method |14| in the complete pf- and Ogg/2-shell. For level density calculations we need to 
solve the interacting shell model problem at finite temperature. The SMMC, based on 
an auxiliary-fields representation of the canonical density matrix at finite temperature, is 
particularly suitable for such calculations. In this letter we discuss only even-even nuclei. 
Odd-A and odd-odd nuclei will be considered in a future publication. For an appropriately 
chosen interaction, even-even nuclei are not subject to the Monte Carlo sign problem, and 
thus their level densities can be calculated more accurately. Except for even-odd effects due 
to the pairing correlations (which are reflected in the backshift parameter A), important 
features of the level density systematics can already be deduced from the study of even-even 
nuclei. 

We use the isoscalar Hamiltonian discussed in Ref. ||. The single-particle energies are 
calculated in a Woods-Saxon potential plus spin-orbit interaction The two-body residual 
interaction includes monopole isovector pairing whose strength is determined from experi- 
mental odd-even mass differences, and a self-consistent surface-peaked interaction |15| . The 
latter is expanded in quadrupole, octupole and hexadecupole terms that are appropriately 
renormalized. The Hamiltonian is uniquely determined for each nucleus by the above proce- 
dure, with no adjustable parameters remaining. Our Hamiltonian includes the dominating 
collective components of realistic effective interactions \T&, yet has the advantage that it 



satisfies the modified sign rule ]T7|, and therefore has a good Monte Carlo sign for even-even 



nuclei. This enables us to perform accurate Monte Carlo calculations down to temperatures 
that are low enough to extract reliable ground state energies. The method used to extract 
the total level density is explained in Ref. ||. The canonical thermal energy {H)p (at a fixed 
proton and neutron number) is calculated in SMMC as a function of inverse temperature 
P, and then integrated to find the canonical partition function Z(f3). The level density p is 
evaluated in the saddle-point approximation 

p(E) = (2nf3- 2 C)- 1/2 e s , (1) 

in terms of the canonical entropy S(E) = (3E + lnZ and the heat capacity C = —(3 2 dE/df3. 
Here the relation between energy and temperature is determined by the saddle-point con- 
dition E(/3) = (H)p. The parity-dependence of the level density is calculated using the 
parity-projection method described in Ref. |§. 

To compare with experimental data, it is necessary to find p as a function of the excitation 
energy. It is thus important to determine accurately the ground-state energy E : any error 
in E will directly affect the value of the backshift parameter. The ground state energy 
can be obtained from the thermal energy in the limit of zero temperature, i.e. E(j3 — > oo). 
Since the auxiliary-fields propagator matrix used in the Monte Carlo method becomes ill- 
conditioned [in| at large (3 (^ 5 MeV -1 ), we have to determine the ground state energy from 
finite (3 calculations, where the lowest excited states still contribute to the thermal energy. 
The ground state of even-even nuclei is J p = + , while the first excited state of almost all 
even-even nuclei is J p = 2 + with an excitation energy of ~ 1 MeV in the 50 ^ A < 70 
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mass region. For sufficiently low temperatures {(3 ^ 2 — 3 MeV -1 ), the main contribution 
to E{(3) — E arises from thermal excitation to the 2\ state. In the two-state model of Of 
and 2f, we have 

E(P)nE + V E x {2l), (2) 

where E x (2f) is the excitation energy of the 2f and r\ = 1/{1 + exp[j3E x (2f)]/5} (where 
the factor 5 accounts for the spin degeneracy of the 2 + state). The ground state energy E 
can then be extracted by a two-parameter fit of (0) to the large-/? Monte Carlo data for 
E(f3). However, in SMMC we can also calculate (J 2 )/3, the canonical expectation value of 
J 2 , where J is the angular-momentum operator. In the two-state model 

(J*)p « 6r? . (3) 

By fitting the Monte Carlo data of {J 2 )p to we can first extract E x (2f), and then use this 
value in (fj) to determine the ground state energy E$ more accurately by a one-parameter 
fit. In Fig. |I|, the extracted values of Eq are shown for 58 Fe as a function of f3. The present 
procedure gives a stable E value beyond p = 2.5 MeV -1 . We adopt the E and E x (2f) 
values by averaging the values in the range (3 = 2.5 to 3 MeV -1 . Our method not only 
provides a more accurate determination of E but also gives E x (2f), although with less 
accuracy (typically ~ 0.1 — 0.2 MeV) than the ground state energy (~ 0.05 MeV). In Fig. |2| 
we compare the calculated excitation energies with the experimental data, and a fairly good 
agreement is observed. This is another confirmation that our present interaction includes 
properly the dominating collective features of the realistic nuclear force. We remark that 
the first observed excited state in 72 Ge is a + , an exception to the usual 2 + . It is not 
clear whether our present Hamiltonian can reproduce this O2" state. Experimentally, the 2f 
state lies only 0.14 MeV above the Of- Because of the spin degeneracy factor, the thermal 
weight of the 2f is still about three times larger than that of 0^ for j3 = 3 MeV -1 . Hence 
we can neglect the contribution of this 0^~ in the above procedure. The reliability of our 
Hamiltonian was further confirmed by the average energy of the mass quadrupole excitation 
in 56 Fe [|. 

The BBF level density [] at excitation energy E x is given by 

P(E X ) « Q^a-\{E X - A)-ieV^-A) (4) 

with g = 2 for the total level density and g — 1 for the parity-projected densities. The 
SMMC total level densities are well described by (f|), and are in good agreement with the 
level densities that are reconstructed from experimentally determined a and A. For example, 
in Fig. we compare the SMMC total level densities of 60 Ni and 68 Zn with the experimental 
ones HH and find good agreement, similar to that found in Ref. for the 56 Fe level 



1 The density given by Eq. (Q) is sometimes referred to in the literature as "state density" where 
each level with spin J is weighted by a factor of 2 J + 1 (to include its magnetic quantum number 
degeneracy). All densities calculated in this paper are state densities. 
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density. The calculated parity-projected level densities are also well fitted to the BBF if 
parity-specific values for a and A (denoted by a± and A±) are used. In general these values 
differ from those for the total level densities. 

We have extracted the level density parameters a and A from SMMC level density 
calculations for a number of /^-stable even-even nuclei in the 50 ^5 A < 70 region: 54_58 Fe, 
58-64^^ 64 - 70 Zn and 70,72 Ge. In some of the previous analyses of level densities (see e.g. 
Ref. 0), a simple empirical function was assumed for the backshift parameter A, and the 
single-particle level density parameter a was then fitted to the experimental data. In Ref. |J 
both a and A are fitted to the experimental data. Empirical formulae for a were proposed in 
Refs. 0)0]) but different formulae yield quite different level densities. Furthermore, there still 
remains a discrepancy between the experimental and the various empirical values of a. While 
these empirical approaches describe the global systematics of the level densities through the 
nuclear periodic table |§, important nuclear structure effects may be overlooked, resulting 
in inaccurate level densities for at least some of the nuclei. Our present calculations seem to 
reproduce available experimental densities with good accuracy (typically within a factor of 
2) for E x < 20 MeV. Consequently, our extracted level density parameters are expected to 
be more accurate than the empirical ones. Our values for a and A are obtained by fitting 
the BBF to the SMMC results in the energy range 4 < E x < 22 MeV. 

Figures |] and |5| show the calculated values of a and A, as well as the parity-dependent 
ones a± and A±, as a function of A. Isotopes are connected by dotted lines. For compar- 
ison, the values obtained from the empirical formula of Ref. for the total level density 
parameters are shown by solid lines. Even-odd staggering in A due to pairing correlations is 
not observed here since only even-even nuclei are considered. It is interesting to determine 
whether shell effects can be observed in the level density parameters. Among the nuclei 
studied, 54 Fe and the Ni isotopes have / 7 / 2 -shell closure for protons (Z = 28) or neutrons 
(N = 28), respectively. In the empirical values of the backshift parameter A of Ref. [[/J, no 
shell effects are seen for these nuclei. In contrast, we find enhancement of A at Z = 28 or 
N = 28, both for the total and parity-projected level densities. Except for the Z = 28 or 
N = 28 nuclei, the present A values for the total level densities are close to those of Ref. [[7| . 
On the other hand, we do not observe shell effects at Z = 28 or N = 28 in the single-particle 
level density parameter a, which increases rather smoothly as a function of A. Compared 
with the empirical estimates of Ref. [|7]], the present values of a are not very different in the 
region A < 65, but the Z-dependence within isobars tends to be weaker. For A > 65 our 
a values are considerably smaller than those of Ref. |7[]. For 68 Zn, for example, we obtain 
a = 7.79 as compared with the empirical value of a = 8.37 0. Our value lies in between 
the experimental values of a = 7.25 (assuming half the rigid-body moment of inertia) and 
a = 7.97 (assuming the rigid-body moment of inertia) |]], and is closer to the rigid-body 
value. We shall return later to the systematics of a. 

We turn next to the parity-dependence of the level densities. For A < 60 nuclei we 
observe that the level density parameters are different for positive- and for negative-parity 
levels; while a + is close to a of the total level density, we find that a_ is larger than a+. 
As A increases, a_ approaches a + (and therefore a). Similarly, for A < 60 nuclei A_ is 
substantially larger than A + and both A± are different from A. The difference between 
A + and A_ becomes smaller as A increases. The observed parity-dependence originates in 
the subshell structure in this mass region, where negative-parity states in even-even nuclei 
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are possible only when the g$/ 2 level is populated. Because of the energy gap between the 
pf and gg/2 orbits we expect the negative-parity level density to be lower than the positive- 
parity level density at low energies. Thus the backshift A_ should be larger than A + . 
On the other hand, at high excitation energies positive- and negative-parity level densities 
are approximately equal, as in the Fermi gas model. Therefore, in the low energy region 
the negative-parity density is expected to rise more quickly as a function of energy , i.e. 
a_ > a + . As A increases, excitations (especially of neutrons) to the gg/2 orbit become easier, 
lessening the difference between p + and p_. For A > 65 we find that the values of a± and 
A± are close to those of the total level densities. 

In the conventional Fermi gas model the single-particle level density parameter a is 
predicted to be proportional to A Q. To investigate the A-dependence of a, it is customary 
to define the quantity K = A/a. In the Fermi gas model K is nearly constant (~ 16 MeV). In 
the empirical model of Ref. [^], K decreases nearly linearly as A increases within a family of 
isotopes, i.e. the level density increases more rapidly as a function of E x for heavier isotopes. 
On the other hand, the empirical formula of Ref. |l[ (see Eq. (10) of [l[]) predicts a gradual 
increase of K. In the present microscopic SMMC calculations, we find that K (for the total 
level densities) depends smoothly on the neutron number N, and is almost independent of 
the proton number Z as is shown in Fig. || K ~ 10 MeV at iV = 28, and decreases towards 
the middle of the N = 28 — 50 shell to a value of ~ 8.5 MeV. This behavior clarifies the shell 
systematics of a. A similar behavior is observed for K + = A/a + . Although we cannot make 
definite conclusions because of the large statistical errors, K_ = A/a- seems to be roughly 
constant. 

In conclusion, using the shell model Monte Carlo method in the complete (pf + Ogg/2)- 
shell, we have calculated microscopically total and parity-projected level densities for (3- 
stable even-even nuclei in the Fe to Ge region. We have studied the systematics and shell 
effects in both the single-particle level density parameter a (as well as the parameter K = 
A /a) and the backshift A. The shell structure observed in A at Z = 28 or N = 28 and the 
systematics of K differ from those given by various empirical formulae. 

This work was supported in part by the DOE grant DE-FG-0291-ER-40608, and by 
the Ministry of Education, Science and Culture of Japan (grant 08740190). Computational 
cycles were provided by the IBM SP2 at JAERI and Fujitsu VPP500 at RIKEN. 



2 Our BBF parametrization of the parity-dependent level densities is valid only for E x < 20 MeV. 
At higher energies p+ ~ p_ and therefore a + ~ a_ should hold. 
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FIG. 1. Thermal energy E(j3) in SMMC (crosses) and ground state energy Eq (diamonds) for 
58 Fe extracted via Eqs. (||) and (||). The solid line is the average of the values for Eq between 
P = 2.5 to 3 MeV" 1 . 
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FIG. 2. Excitation energies E x {2+) of the first 2+ state: comparison between the SMMC values 
(calculated from Eq. (|3|)) (crosses) and the experimental values (diamonds). Mass numbers are 
shown for several nuclei below the symbols. 
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FIG. 3. Total level densities: comparison between SMMC (solid squares with error bars) and 
experimentally determined level densities (solid lines), a) 60 Ni. The experimental BBF parameters 
are taken from Ref. ||. b) 68 Zn. The experimental BBF parameters are taken from Ref. |Q] 
assuming rigid-body moment-of-inertia. 
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FIG. 4. Single-particle level density parameter a as a function of A, obtained from fitting the 
SMMC level densities to the BBF (Q). Isotopes are connected by dotted lines. Left: a for the 
total level densities (squares) in comparison with the values obtained from the empirical formula 
of Ref. (solid lines). Right: a+ (circles) and a_ (triangles) of the positive- and negative-parity 
SMMC level densities, respectively. 
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FIG. 5. Backshift parameter A as a function of A, obtained from SMMC calculations for total 
(left) and parity-projected (right) level densities. Conventions are as in Fig. ||. 
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FIG. 6. -fT = A/a as a function of neutron number N from the SMMC calculations for total 
(left) and parity-projected (right) level densities. Isotopes are connected by dotted lines. 
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